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We carry out Faddeev calculations of three-alpha (3a) and two-alpha plus A (aaA) systems, using 
two-cluster resonating-group method kernels. The input includes an effective two-nucleon force for 
the act resonating-group method and a new effective AN force for the Aa interaction. The latter 
force is a simple two-range Gaussian potential for each spin-singlet and triplet state, generated 
from the phase-shift behavior of the quark-model hyperon-nucleon interaction, fss2, by using an 
inversion method based on supersymmetric quantum mechanics. Owing to the exact treatment of 
the Pauli-forbidden states between two a-clusters, the present three-cluster Faddeev formalism can 
describe the mutually related, aa, 3a and aaA systems, in terms of a unique set of the baryon- 
baryon interactions. For the three-range Minnesota force which describes the aa phase shifts quite 
accurately, the ground-state and excitation energies of ^Be are reproduced within 100 ~ 200 keV 
accuracy. 

PACS numbers: 21.45.+V, 21.30.-x, 13.75.Cs, 12.39.Jh 



I. INTRODUCTION 



In spite of much effort to incorporate microscopic fea- 
tures of the alpha-alpha (aa) interaction, a consistent 
description of the three-alpha (3a) and two-alpha plus 
A (aaA) systems has not yet been obtained so far in 
the Faddeev formalism. The most favorable description 
of the aa system is the aa resonating-group method 
(RGM) [1]. Although some of the previous works deal 
with the aa RGM kernel explicitly in the 3a-cluster 
Faddeev formalism, they usually yield a large overbad- 
ing for the ground state and sometimes involve spurious 
states because of an incomplete treatment of the Pauli- 
forbidden states in the 3a model space [2-5]. Various 
types of 3a orthogonality condition models (OCM) [6- 
8] also yield a similar overbinding for the ground state, 
although the effect of the Pauli principle between a clus- 
ters is satisfactorily treated in each framework. Only one 
exception for this rule is the 3a OCM in Refs. [9, 10], in 
which the Pauli forbidden components described by the 
aa bound-state solutions of the deep Buck, Friedrich, 
and Wheatlcy (BFW) potential [11] is completely elimi- 
nated. The result is rather similar to the traditional 3a 
Faddeev calculation using Ali-Bodmer phenomenological 
aa potential with a repulsive core [12]. In these cal- 
culations, the ground-state energy of the 3a system is 
less than 1.5 MeV, and a simultaneous description of the 
compact shell-model like ground state and the excited + 
state with well-developed cluster structure is not possi- 
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ble. The origin of the different conclusions in Refs. [8] 
and [9, 10], is spelled out in Ref. [13], in which the exis- 
tence of almost forbidden Faddeev components inherent 
to this 3a OCM using the bound-state Pauli-forbidden 
states of the BFW potential is essential. 

A possible resolution of this overbinding problem of 
the 3a model is found in our new three-cluster Faddeev 
formalism, which uses singularity-free T-matrices (RGM 
T-matrices) generated from the two-cluster RGM kernels 
[14]. In this formalism, solving the Faddeev equation 
automatically guarantees the elimination of the three- 
cluster redundant components from the total wave func- 
tion. The explicit energy dependence inherent in the ex- 
change RGM kernel is self-consistently treated. We first 
applied this formalism to the three-dineutron and 3a sys- 
tems, and obtained complete agreement between the Fad- 
deev calculations and variational calculations using the 
translationally invariant harmonic-oscillator (h.o.) ba- 
sis [14, 15]. Next, this formalism was applied to a Fad- 
deev calculation of the three-nucleon bound state [16], 
which employs complete off-shell T-matrices derived from 
the non-local and energy-dependent RGM kernels of the 
quark-model NN interactions, FSS [17] and fss2 [18]. 
The fss2 model yields a triton binding energy B t = 8.519 
MeV in the 50 channel calculation, when the np interac- 
tion is employed for all the NN pairs in the isospin basis 
[19]. The effect of the charge dependence of the two-body 
NN interaction is estimated to be —0.19 MeV for the tri- 
ton binding energy [20]. This implies that our result is 
not overbinding in comparison with the empirical value, 
B t cxp = 8.482 MeV. If we attribute the difference, 0.15 
MeV, to the effect of the three-nucleon force, it is by far 
smaller than the generally accepted values, 0.5 ~ 1 MeV 
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[21], predicted by many Faddeev calculations employing 
modern realistic meson-theoretical NN interactions. We 
have further applied this three-cluster Faddeev formalism 
to the hypertriton system [22] , in which the quark- model 
hyperon-nucleon (YN) interactions of fss2 yield a reason- 
able result of the hypertriton properties similar to the 
Nijmegen soft-core potential NSC89 [23]. Most math- 
ematical details for the Faddeev equation, employed in 
this calculation, are given in the present paper. 

Here we apply the present three-cluster Faddeev for- 
malism to the aaA model for 9 Be. This hypernucleus 
plays an important role to study the A7V interaction in 
the p-shell A-hypernuclei. From the early time of the 
hypcrnuclcar study, ^Be is considered to be a proto- 
type of a-cluster structure, in which the two a clus- 
ters form a loosely bound subsystem by the effect of 
the extra A hyperon [24]. Since the YN interaction 
is usually weaker than the NN interaction, this sys- 
tem is suitable for studying a subtle structure change 
of the two-a system from 8 Be. In fact, in addition to 
the l/2 + ground state [25-28] with the A-separation en- 
ergy B A ( A Be) = 6.71 ± 0.04 McV [29], the recent 7 -ray 
spectroscopy [30, 31] has revealed the existence of two 
narrow resonances in the excited states, which are sup- 
posed to be 5/2 + and 3/2 + states generated from the 
small spin-orbit splitting in the weak coupling picture of 
8 Be(2+) x A (spin S = 1/2). From a theoretical point 
of view, this is the simplest non-trivial system which re- 
quires the Faddeev formalism with two identical particles, 
involving three Pauli-forbidden states between two clus- 
ters. Several model calculations were already done with 
various frameworks and two-body potentials. Hiyama et 
al. [7] used the OCM for the aa, 3a and aaA systems 
and discussed not only the ground state of 9 Be, but also 
the spin-orbit splitting of the 5/2 + and 3/2 + states [32]. 
They employed simple three-range Gaussian potentials 
for the A7V interaction [33] based on G-matrix calcu- 
lations of various Nijmegen and Jiilich YN one-boson- 
cxchange-potential (OBEP) models. The Aa potentials 
are generated from these AN effective potentials by the 
folding procedure with respect to the (0s) 4 h.o. wave 
function of the a cluster. They introduced a three-a force 
and adjusted the YN parameters to reproduce the bind- 
ing energies of the 12 C and 9 Be ground states. Filikhin 
and Gal [34] used the Faddeev and Faddeev- Yakubovsky 
formalisms to calculate the A Be and A ° A Be ground states. 
They used the Ali-Bodmer aa potential [12] and the so- 
called Isle potential [35] for the Aa interaction. They 
included only S'-wave in the calculation, and reproduced 
the ^Be ground-state energy correctly. However, if one 
includes higher partial waves the Ali-Bodmer aa poten- 
tial yields overbinding for ^Be by more than 0.5 MeV. 
Oryu et al. [36] carried out a aaA Faddeev calculation 
by using the aa RGM kernel and various types of Aa 
potentials in the separable expansion method. Their en- 
ergy spectrum of ^Be is reasonable, but the treatment 
of the two-a Pauli principle in the aaA system is only 
approximate. Since they neglected the Coulomb force, a 



detailed comparison between their calculated results and 
experiment is not possible. Cravo, Fonseca, and Koike 
[37] performed aaA and aan Faddeev calculations by us- 
ing many aa and Aa potentials with the Coulomb force 
included between the two a particles. From the com- 
parison of the results for the 9 Be and 9 Be systems, they 
found an interesting sign change of the quadrupole mo- 
ments and the magnetic moments for some excited states. 
They also pointed out a possibility of negative parity res- 
onances with ^He + a cluster structure in this threshold 
region. 

Our purpose for the aaA Faddeev calculations using 
aa RGM kernels is threefold. First, we develop a gen- 
eral three-cluster Faddeev formalism with two identical 
clusters, in order to apply it to more complex three- 
cluster systems like the hypertriton interacting via quark- 
model baryon-baryon interactions. In the hypertriton 
system, we have to deal with the ANN-HNN coupled- 
channel system which involves a Pauli-forbidden state at 
the quark level in the AN-T,N subsystem. Since the 
baryon-baryon interactions in the quark model are for- 
mulated in the two-cluster RGM formalism, the present 
three-cluster formalism is most appropriate to correlate 
the baryon-baryon interactions with the structure of few- 
baryon systems. The second purpose is to make a con- 
sistent description of the aa, 3a and aaA systems using 
effective NN and AN interactions. This attempt is be- 
yond the scope of the usual OCM framework and the Fad- 
deev formalism assuming only inter-cluster potentials. A 
comparison of the present 3a results with the fully mi- 
croscopic 3a RGM or GCM [38-40] is useful to examine 
the approximations involved in the present three-cluster 
formalism. The third purpose is to present a tractable ef- 
fective AN force for cluster calculations of various p-shell 
hypernuclei, which is not purely phenomenological but 
derived microscopically from quark-model baryon-baryon 
interactions. In particular, this effective AA^ force should 
be able to reproduce the correct A-separation energy of 
^Hc; B A ( A tte) = 3.12±0.02 MeV. Such a Aa interaction 
is indispensable for, e.g., a AAa Faddeev calculation us- 
ing the quark-model AA interaction [41]. In this paper, 
we derive an effective AN force of two-range Gaussian 
form from the phase-shift behavior of the quark-model 
YN interaction, fss2, by using an inversion method based 
on supersymmetric quantum mechanics [42]. 

This paper is organized as follows. In the next section, 
the three-cluster Faddeev formalism with two identical 
clusters is given, together with expressions to calculate 
the expectation values of the two-cluster Hamiltonian 
with respect to the solutions resulting from the Faddeev 
equation. The procedure to calculate the Aa and aa 
T-matrices is also discussed, as well as the treatment of 
the cut-off Coulomb force employed in this paper. In the 
third section, we first briefly discuss the results of the 3a 
Faddeev calculation, and then those of the aaA Faddeev 
calculation. The final section is devoted to a summary. 
Appendix A gives a brief comment on the rearrangement 
factors of three-body systems with two identical parti- 
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cles. The most general case with explicit spin-isospin 
degrees of freedom is discussed. In Appendix B, we de- 
rive a compact formula to calculate the Ka Born kernel 
for arbitrary types of AN interactions. Energies arc in 
MeV and lengths in fm throughout. 



II. FORMULATION 



A. Faddeev equation for systems with two identical 
clusters 



In order to formulate the Faddeev equation for systems 
with two identical particles, we follow the notation of 
Refs. [43, 44] as much as possible. The Jacobi-coordinatc 
vectors are specified by the permutation (a/37), which is a 
cyclic permutation of (123). For example, the momentum 
vectors for the coordinate system 7 in the unit of h are 
defined by 



mpk a — m a kp 



m a + mp 



1 



^ = M [{ma 
P = k a + kp 



- mp) fe 7 



m 7 (k a + kp)] , 



(1) 



where k a , kp, k 1 are single particle momenta of particles 
a, (3, 7 with the masses m a , mp, m 7 , respectively, and 
M = ma+mp+m^ is the total mass. Three different sets 
of the Jacobi coordinates, (p 1 ,q 1 ), (P2,Q 2 ), an d (^3,^3)7 
are related to each other in the standard relationship for 
the rearrangement. We choose the coordinate system 7 = 
3 as the standard set of Jacobi coordinates and assume 
that particles 1 and 2 are the two identical particles with 
a common mass mi =ni2. We incorporate the symmetry 
property for the exchange of particles 1 and 2 into the 
Faddeev formalism by assuming the total wave function 
*(p,g) as 



*(P 3 > I3) = V>(?>3> Is) ± V(-Pl) <?l) + </>(/>2> 92) 
with lp(-p 3 , <? 3 ) = ±V>(P3; 9s) i 



(2) 



where the upper (lower) sign is applied for identi- 
cal bosons (fermions). The requirement p, q) — 
±*(p, q) is satisfied from this ansaz. 

In the application to the aaA system, two a clusters 
are numbered 1 and 2, and the A hyperon is numbered 
3. Since the technique to handle the rearrangement of 
the Jacobi coordinates in the Faddeev formalism is well 
known [43] , we only give the specification scheme of chan- 
nels and the final Faddeev equation after partial-wave 
decomposition. We give expressions both in the LS- 
coupling and jj-coupling schemes for later convenience. 
For the ^Be system, the 7 channel is specified by 7 = 3 
with (jaf3) — (312) in Eq. (1). A set of quantum num- 
bers in the 7-channel is specified by 7 = [(X£)L1/2]JJ Z 
in the LS'-coupling scheme and [X(£l/2)j] JJ Z in the jj- 
coupling scheme with the angular-spin wave functions 



(L5-coupling) 



[YxipMiq,) £i(3)kU ^-coupling) . 



(3) 



Here, Y {xe)LLz (p,q) = [Yx(p)Y e (q)] LLz , £i(3) is the spin 
wave function of A, and <j) a is the internal wave function 
of the a particle. Similarly, we define the (3 channel by 
(3 = 2 with {(3^a) = (231), and a set of quantum numbers 
(3 = [(£ 1 £ 2 )Ll/2}JJ z (LS-coupling) and [(£ 1 1/2)I£ 2 \JJ Z 
(jj-coupling) with 



[ [>w 2 )i,(P2><y £i(3)]j./, 



(LS'-coupling) 



[Y1AP2) U{3)]iY ea (q 2 )]j Jt (jj-coupling) . 



(4) 



The a channel is specified by a — 1 with (a(3j) = (123), 
and the quantum numbers similar to those of the (3 chan- 
nel. The partial-wave decomposed Faddeev equation for 
the two components tp and (p in Eq. (2) is given by 
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4M 



A' 
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q' Z dq' / dx (p\T x E - 



ft 2 8 + C 
AM N 4( 



x J2^x9^(q,q;x)-^ FI ipp(p 2 ,q) , 



^ 
vp(p, q) 



(5a) 
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8M N V C 4 + C 
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h 2 8 + C 2 
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8M N 4 + C 
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(5b) 
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where ( = (M^/Mm) is the mass ratio of A to the nucleon and 



Pl = P (q',y-,x) ip' 1 =p(q', 1 ^q;x) 
P2=p{q,£- C q';x) ' \p' 2 =p(q,y';x) 









I Pi 


= p 


q', 4Tc<?;^j 








[ Pi 


= p 





(5c) 



with p{q,q']x) = \J q 2 + q' 2 + 2qq'x. The T-matrices, T\ and Tg 1 , are discussed in Subsecs. II. D and II. C. The 
rearrangement factors for the ip - <p or ip - -0 cross terms are given by 



g-yp(q,q';x) = gp 7 (q',q;x) = q x '^ +X2 q 

Ai+A 2 = A Ai + A^=^i 



,Ai+A 2 ( 1 

2 



(irc) A2 E( 2fc + 1 )4f lfe ^W' (e) 



where Pk{x) is the Legendre polynomial of rank k. The reduced rearrangement factor g*p Xlk is expressed as 

(LS— coupling) 



Ai X[k 



E L (-V I+J+L+ei+1 3T(Ly 
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i 




{ J L \ \ 




j 







G 



Ai \' r kL 



(ij-coupling) 



(7) 



with j = \/2] + 1 etc., and the spatial angular- momentum factor ^(aV) 1 ^^) ' n Appendix Eq. (A9). For the ip - if 
type rearrangement, these factors are given by 



9pp,(q,q';x) = £ ]T ^i+V^ 

Ai+A 2 =^i Ai+A 2 =«i 



4 + C 



A 2 +A 2 



^(2fc+l).g^ fe P fe ( 



(8) 



with 



Ai At k 



J L \ \{ J L \ 
i 2 \\ £[ I' i' 2 



G 



Ai A^ kL 

(« 2 ),(^ 2 ) 



(LS— coupling) 
{jj- coupling) 



(9) 



B. Calculation of e 7 and ep 

In this subsection, we derive some formulas to calculate 
expectation values of the two-cluster Hamiltonians, /i 7 + 
V^j RGM (e 7 ) and hp + Vp, where h 7 is the kinetic-energy 
operator of the 7-pair etc. In the present application, 
K Y RGM (e 7 ) is the aa RGM kernel and Vp is the Aa kernel. 
We deal with the energy dependence of the aa RGM 
kernel self-consistently by calculating 



(*|^ 7 + K KGM (e 7 )|*) 



(10) 



for the normalized Faddeev solution <!/. The potential 
term of the matrix element in Eq. (10) is most easily 
obtained from various matrix elements of the kinetic- 
energy operators. Suppose $ be a sum of three Fad- 
deev components, ^ = ip a + ipp + ipj. Then the Faddeev 
equation (E - H )ip 7 = with V" 7 = V 7 RGM (£ 7 ) and 
H = h 7 + /i 7 yields (4'|V r 7 |4') = (ip 7 \E - H \V). Thus 
Eq. (10) becomes 

e 7 = £7<^ 7 |*> - (V 7 |ff |*> + (*|/fcrl*) • (11) 

We can write a similar equation also for the (3 pair. We 
calculate ep, although the self-consistent procedure is not 



necessary for the Aa interaction. The kinetic energy term 
{^\hp\^} is obtained from (*|ft 7 |*) as follows. Using the 
momentum Jacobi coordinates in Eq. (1), we can easily 
show 

(rap + m 7 )h a + (m 7 + m a )hp + (m a + mp)hj 
= MH . (12) 

For two identical particles with m a = rap, this relation- 
ship yields 

h a + hp = ^ H — h-/ , (13) 

rap + m 7 rap + m 7 



and 



(*IM*> 



M 



2(mp + m 7 ) 
rap 



rap + m 7 

Thus we find, for the aaA system, 



<*|H |*> 

^1^1*) . (14) 



ep = Etyp\9) + ^±^fa\H \*) 

[<^|ff |*> - <*N*>] ■ (15) 



4 + C 
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We need to calculate the overlap matrix elements (-0 7 1 (^\h^\ip^} + 4 (V> 7 |/i 7 |V'/j) + 2 ^h^tf) p + ip a ). These 
{tppYi), and (tpj\HQ\^>) = (ip 7 \H \tp 7 ) + 2 (^> 7 |if |V'/j)) are calculated from ?/> and <p by using the recoupling tech- 
(^g|.ffo|*) = {■f/'-y |i^o|"0/3} + ("0/3 l^o|V , /3 + , W^I'I') = niques developed in Appendix A. The final result is 



f°° h 2 ( 

= E/ p 2d p q2dq m^ [p 2 

+ T2 q 2 dqq' 2 dq' dx -0 7 (pi, g)— — ( pi 2 + — — <? 2 ) —rg^piq, q'l x) — j-tpfj(p2 , q') , (16a) 

(^iff i*> =n[p 2d p^ 2d ^ {^r p2 + Trc q2 ) [ ^ q)]2 

1 /*°° 2J /2, , r 1 , ,~ v ^ /4 + c~2 8 + c 2 \ i . , . i 

+ 2^J q qq q J 1 Vfi ^' q >8M^\C~ Pl + 4 + C 9 J ~7T^/3'(9.9 ; ^-377 </W2, 9 ) 

l . z* 00 2 f 1 h 2 / 8 + C \ l l 

+ oE/ q 2 dqq' dq' dx^ J (p 1 ,q)—^(p 1 2 + ——q 2 )^g 1 p(q,q';x)^-(p f 3(p 2 ,q'). (16b) 
The overlap integrals are obtained by setting i? — > 1- Furthermore, (\l/|/i 7 |\I/) is given by 

/■oo t2 

<*N*> = E/ o p 2dpq2dq m^ p2 

roo 9 f 1 H 2 1 1 

+ E / 9 2 rf<?<?' dq' / dxtP 1 (p 1 ,q)——p 1 2 ^-g lf 3(q,q';x)— r -(p (p2,q') 
J-l 2M N Pl A P2 l1 

+ £i ^ P ^^' 9) 8^{[ P2+ (4T|) 

+ E / q 2 dqq' 2 dq' dxip^p^q) (q 2 + q' 2 -2qq'x) —j-g^^q, g ' ; x )— J - 7 ip l3 ,(p 2 , q') . (17) 

J^Jo J-i lbM N V / pi 1 pj'i 

Here, is given by 

//J/J' = y dp 2 dq 2 (/3\p 2 q 2 )(p 2 -q 2 )(p 2 q 2 \P') = (-1)^+^F 1 ^(^ 1 010|£ 1 ') (^OlOl^') 

(-1) L { £/ i } (LS-coupling) 



EJ-D— /T'(£) 2 { I \\ } { ^ L p | } { £, £ f } (ii-coupling) 



(18) 



I 

C. Aa T-matrix and effective AN potentials where fj, — [4£/(4 + ()]Mn is the Aa reduced mass and 

E = — (/i 2 /2/i)7 2 is a negative energy. The partial- wave 
components V((p,p') for the Aa Born kernel V(p,p') are 
The Aa T-matrices are obtained by solving the defined through 
Lippmann-Schwinger equation 

V(p, P ') = Vt(p,p') J2 Y tm(pyYirn(p') , (20) 

4tt 2^ f°° em 
T e (p,p>; E) = V t M) j^ Y2 I k 2 dk V t (p, k) ^ in Eq (5b) ig related tQ w . E) 

1 , with an extra factor 47r/(27r) 3 . 

x 7 2 + k 2 Tl ( fc ' P ' E > ' ( 19 ) For the effective AN potential, we assume a Minnesota- 
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type central force [45] 



VAN 



vCE) 



I -Pa 



v{ 3 E) 



1 + Pa 



u 2 — u „ 

- H P r 

2 2 r 



(21) 



where u( 1 i?) and v( 3 E) are simple two-range Gaussian 
potentials generated from the 1 5o and 3 Si phase shifts 
predicted by the quark-model A7V interaction, fss2. We 
use the inversion method based on supersymmetric quan- 
tum mechanics, developed in Ref. [42], to derive phase- 
shift equivalent local potentials. These potentials are 
then fitted by two-range Gaussian functions. These are 
given by 



vCS ) 



-128.0 cxp(-0.8908 r 2 ) 
+1015 cxp(-5.383 r 2 ) , 
-56.31/ cxp(-0.7517 r 2 ) 
+1072 cxp(-13.74 r 2 ) , 



(22) 



where / = 1 and r is the relative distance between A and 
N . In the following, we call this effective AN potential 
the SB potential. Figure 1 shows that these potentials 
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FIG. 1: AN-T.N (a) and 3 Si (b) phase shifts for the 
isospin 1—1/2 channel, calculated with fss2 [46] (solid and 
dashed curves) and with the SB potential (circles). 
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TABLE I: Aa spin-flavor coefficients for the Minnesota-type 
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fit the low-energy behavior of the 1 S'o and 3 Si AA^ phase 
shifts obtained by the full AN-Y.N coupled-channel RGM 
calculations of fss2. In the 3 Si state, only the low-energy 
region is fitted, since the cusp region cannot be fitted in a 
single-channel calculation. This potential overestimates 
the 3 Si phase shift in the higher energy region. The 
procedure to calculate the Aa Born kernel for the simple 
(0s) 4 a-cluster wave function is discussed in Appendix 
B. Here we only give the final result for the partial- wave 
components: 

4 

Vt{q f , Qi) - J2 i X d V i ft! K i) + X e V e (if' If, Ki)] . 



(23) 



Here, X l d and X\ are spin-isospin factors defined in Eq. 
(B10) and tabulated in Table I for the present two-range 
Gaussian potentials. The explicit functional form of 
Vf {q/, qi] Ki) and V[ (q f , q t ; Ki) are given in Eq. (Bll). 

In this paper, we also examine the AN effective forces 
[33] used by Hiyama et al. [7] for comparison. These 
potentials are generated from the G-matrix calculations 
of various OBEP potentials. They are parameterized as 



VAN 



i=l 



,(0 



,(*) 



even 1 aa even 



(<Tl • <T 2 ) 



+ 



1 1 ) (i) 
V odd + V av odd (o-l -0-2) 



1 - P r 



l+Pr 

2 

-('•/ft) 2 



(24) 



Since the spin-spin term does not contribute to the spin 
saturated a-cluster, the spin-isospin factors in Eq. (23) 
(with Ki — ► l/(/3i) 2 , i = 1 - 3) are given by 



A d — L [Vo oven + V odd J 

X'' - 2 /V J) -7; (i) ^ 

A e — z \ u oven u odd J 



(25) 



The explicit values for d^ cvc11 and v^' odd (i = 1 - 3) gen- 
erated from Nijmcgcn models, NS, ND, NF, and Jiilich 
potentials, JA, JB, are given in Ref. [7]. [Tabic V of 
Rcf. [7] includes a misprint for NS: the width parameters 
(3i for this potential are 1.50 - 1.0 - 0.55, instead of 1.50 
- 0.90 - 0.50 for the other potentials.] 

The binding energy of the Aa bound state depends 
on the h.o. width parameter v of the en-cluster. Table 



7 



TABLE II: Bound-state energies for the Aa system, E(\¥Le) 
(in MeV), calculated by the original SB potential with / = 
1. The h.o. width parameters, v = 0.275 fm~ 2 and v = 
0.257 fm~ 2 are assumed for the (0s) 4 a-cluster. The experi- 
mental value is £ oxp (iHe) = -3.12 ± 0.02 MeV. 



u 


v = 0.275 fin" 2 


v = 0.257 fm~ 2 


1 


-4.975 


-4.747 


0.6 


-4.946 


-4.728 



II shows that the SB potential of Eq. (22) overbinds the 
^He energy by more than 1.6 MeV. It also shows that the 
u-dependence is very weak, which implies that ^He is an 
S-wave dominated system. It is well known that a central 
single-channel AN effective force that fits the low-energy 
AA total cross section and the ground-state energies of 
^H, and ^He always overestimates the ^He binding 
energy by more than 2 MeV, due to a lack of A-E mix- 
ing and the tensor force [47-50]. In order to circumvent 
this difficulty, we introduce a reduction factor / in the 
attractive part of the 3 Si potential in Eq. (22) for the fol- 
lowing Faddeev calculations. The choices / = 0.8821 for 
v = 0.275 fm~ 2 and / = 0.8923 for v = 0.257 far 2 re- 
produce the desired value S(^He) = —3.120 MeV, when 
the pure Serber type SB potential with u = 1 is used. 
The Aa bound-state energies predicted by the NS - JB 
effective AA potentials deviate from the original fit in 
Ref. [7] by 110 keV - 170 keV (-3.23 - -3.29 MeV). 
This is because they used a slightly different expression 
from ours for the exchange term of the Aa potential. For 
the aaA Faddeev calculations using the Minnesota three- 
range force for the aa RGM kernel (y = 0.257 far 2 ), we 
readjusted the strength of the original NS - JB AA?" po- 
tentials in order to fit the precise ^He energy, —3.120 
MeV. This is achieved by slightly (less than 0.36%) mod- 
ifying the strength of the short-range repulsive term (the 
third component) of the original G-matrix potentials. 

The Aa phase shifts are also calculated, although there 
is no cxpcrimetal information. The S'-wave phase shift 
shows a monotonic decrease from 180°, similar to Fig. 
9 of Ref. [36]. In the energy region T c . m .(Aa) = 0-20 
MeV, the phase shifts of the higher partial waves rapidly 
decrease, starting from 20° - 30° for the P wave. This 
implies that the Aa potential is very much of the Wigncr 
type, and our lack of knowledge of the Aa interaction in 
higher partial waves may not become a serious problem 
in the Faddeev calculations. 



D. oca T-matrix and effective NN Potentials 

The aa T-matrices used for the 3a and aaA Fad- 
deev calculations are generated from the aa RGM kernel 
which uses an effective AW potential similar to Eq. (21). 
In the notation used in Ref. [51], the aa RGM kernel, 



V RGM (e) = Vd + r D CL + G + G CL +sK consists of the 
direct potential Vb, the direct Coulomb potential Vd° L , 
the sum of the exchange kinetic-energy and interaction 
kernels, G = G K + G v , the Coulomb exchange kernel 
G CL , and the exchange normalization kernel K . We have 
to eliminate redundant components from the energy- 
dependent partial wave T-matrices, T\(p,p'; E, e), which 
satisfy the Lippmann-Schwingcr equation similar to Eq. 
(19). This is necessary only for the S'-wave (A = 0) and 
.D-wave (A = 2) components, for which there exist two 
and one h.o. Pauli redundant states, u nX (p), respectively. 
Here, u n \(p) are essentially the h.o. wave functions in 
the momentum representation with the total h.o. quanta 
N = 2n + A = and 2, satisfying Ku nX = u nX . They are 
explicitly given by 



(2tt) 3 / 2 

u nX {p) = {-l) nK —±=- Rnx[p, . , 

V47T V / 



(26) 



with [i = 2, in terms of the standard three-dimensional 
h.o. radial wave function R n g(r, v) in the coordinate rep- 
resentation. The RGM T-matrices defined in Ref. [14] 
are calculated by 



T x (p,p';E,s) = T x (p,p';E,e) 

h 2 ( 7 2 +p 2 ) (l 2 +p' 



+ 



AM N (7 2 + k 2 ) 
for A=/° , 



El=o u «o(p) u n0 (p') 
(27) 



where 7 2 = -(4M N E/h 2 ) and k 2 = {AM N e/h 2 ). 
For the higher partial waves with A > 4, we define 
f x (p,p';E,e) = T x (p,p';E,s). The RGM T-matrices in 
Eq. (27) satisfy the orthogonality condition 



u n \(p) 

4-7T 4M N 
= (2tt) 3 h 2 



/ p' 2 dp'T x (p,p';E,e) 
Jo 



u nX (p') 
l 2 +p' 2 



(28) 



for n = 0, 1 (A = 0) and n = (A = 2). Owing to 
this relationship, we can prove the orthogonality of the 
total wave function Eq. (2) to the Pauli-forbiddcn states 
u nX {p). 

For the effective AW force, we mainly use the three- 
range Minnesota (MN) force [45] with the exchange- 
mixture parameter, u = 0.94687, and the h.o. width pa- 
rameter, v — 0.257 far 2 , for the (0s) 4 a-clusters. We 
also use the two-range Volkov No. 1 (VN1) and No. 2 
(VN2) forces [52], in order to compare our 3a results 
with the microscopic RGM [38, 53] and GCM [39] cal- 
culations. The Majorana parameters m of the Volkov 
forces and the h.o. width parameters are m = 0.575 
and v = 0.2515 far 2 for VN1, and m = 0.59 and 
v = 0.275 far 2 for VN2. The aa RGM calculations us- 
ing these effective AW forces and the complete Coulomb 
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kernel reasonably reproduce the empirical aa phase shifts 
of the 5-, D-, and G-waves, as well as the .S'-wave reso- 
nance near the aa threshold. However, the best fit to the 
experiment is obtained by the three-range MN force. For 
the VN2 force, the .S'-wave resonance appears as a bound 
state with the binding energy B aa = 245 keV. Although 
the VN1 force reproduces this resonance, the overall fit 
to the aa phase shifts is less impressive compared to the 
MN force. In the RGM calculation, the precise determi- 
nation of the resonance energy is not easy even in the 
two-a system, because of the presence of the Coulomb 
force. In the present Lippmann-Schwinger formalism in 
the momentum representation, the method by Vincent 
and Phatak [54] is used for solving the scattering prob- 
lem with the full Coulomb force at the nucleon level. We 
find that the + resonance energy is 0.18 MeV and 0.14 
MeV for VN1 force and the MN force, respectively. This 
should be compared with the experimental value 0.092 
MeV. 

For the Coulomb force in the 3a and aaA Faddeev cal- 
culations, we use the cut-off Coulomb force at the nucleon 
level 

^(r) = ^^ f^-r), (29) 

with the cut-off radius Rc, although an exact treat- 
ment of the point Coulomb force exists for bound-state 
nuclear three-body problems with two charged particles 
[55]. Here 8(x) is the Heaviside step function. For the 
most compact 3a ground state, this approximation with 
Rc = 10 fm is good enough to obtain 1-2 keV ac- 
curacy. The Coulomb exchange kernel for Eq. (29) is 
calculated analytically. The partial-wave decomposition 
of the aa RGM kernel is carried out numerically using 
the Gauss-Legendre 20-point quadrature formula, when 
the Coulomb force is not included. When the cut-off 
Coulomb force with Rc = 14 fm is employed, it is in- 
creased to the 30-point quadrature formula to obtain an 
accuracy within 1 keV for the exchange Coulomb kernel. 
The direct Coulomb term is separately integrated with 
enough number of numerical integration points. 

III. RESULTS 



TABLE III: Results of 3a Faddeev calculations, using the 
aa RGM kernel, with and without the Coulomb effect. The 
parenthesized numbers indicate the results when the cut-off 
Coulomb force with Rc = 10 fm are included at the nucleon 
level. Partial waves up to A max are included in aa and (2a)- 
a channels. The heading £20, is the expectation value of the 
two-a Hamiltonian with respect to the 3a bound state solu- 
tion, E;i a the 3a bound-state energy, and C( 4) the overlap 
between the 3a bound-state wave function and the SU3 (04) 
shell-model configuration. For the MN force, the result of 
the variational calculation using the translationally invariant 
h.o. basis (h.o. var.) is also given for comparison, where h.o. 
quanta up to N = 60 are included. 



Force Amax £2a ^(04) 

4 9.657 (10.887) -10.751 (-5.206) 0.900 (0.879) 

VN1 6 9.531 (10.779) -10.926 (-5.365) 0.896 (0.875) 

8 9.530 (10.778) -10.927 (-5.366) 0.896 (0.875) 

4 8.583 (9.608) -11.202 (-5.781) 0.826 (0.795) 

VN2 6 8.449 (9.505) -11.415 (-5.967) 0.821 (0.790) 

8 8.447 (9.503) -11.417 (-5.969) 0.821 (0.790) 

4 12.032 (13.603) -15.616 (-9.433) 0.979 (0.973) 

MN 6 11.905 (13.482) -15.777 (-9.591) 0.978 (0.971) 

8 11.904 (13.481) -15.779 (-9.592) 0.978 (0.971) 

h.o. var. 11.903 (13.480) -15.781 (-9.594) 0.978 (0.971) 



tion. These 713 points for p > 6 fm - are, however, not 
included for solving the Faddeev equation, since it causes 
a numerical instability for the interpolation. The mo- 
mentum region q = 6 fm -1 - 00 is also discretized by 
the 7J3 point formula just as in the p discretization. We 
take ni-n2-ri3=15-10-5 for the 3a system and 10-10-5 for 
the aaA system, respectively, unless otherwise specified. 
The modified spline interpolation technique developed in 
[56] is employed to generate the rearrangement matrices. 
For the large-scale diagonalization of non-symmetric ma- 
trices, the Arnoldi-Lanczos algorithm developed in the 
ARPACK subroutine package [57] is very useful. 



To solve the Faddeev equation, we discretize the con- 
tinuous momentum variable p (q) for the Jacobi coor- 
dinate vectors, using the Gauss-Legendre m-point (n 2 - 
point) quadrature formula, for each of the three inter- 
vals of - 1 fm -1 , 1 - 3 fm" 1 and 3 - 6 far 1 . The 
small contribution from the intermediate integral over 
p beyond p = 6 far 1 in the aa T- matrix calculation 
is also taken into account by using the Gauss-Legendre 
ri3-point quadrature formula through the mapping p = 
Po + tan{7r(l +x)/4}. We need n\ > 10 and n 3 = 5, so 
that 35 points are at least necessary to follow up the inner 
oscillations of the two-a bound-state wave function and 
the necessary T-matrices for solving the Faddeev equa- 



A. 3a Faddeev calculation 

In order to make sure that our Faddeev equation 
is solved correctly, we first carried out the standard 
3a-particle Faddeev calculation by using the angular- 
momentum dependent Ali-Bodmer potential of d type 
(ABd). We find that the 3a energy, E 3a = -6.423 MeV 
without Coulomb force, is consistent with previous cal- 
culations [9]. Here, we used h 2 /M a = 10.4465 MeV • fm 2 
and e 2 = 1.44 MeV • fm for comparison. When the cut- 
off Coulomb force is included, our value —1.527 MeV 
is 4 keV lower than the —1.523 MeV given in Table 1 
of Ref. [9]. This difference is due to a slightly different 
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treatment of the Coulomb force between the two calcula- 
tions. The small 3a binding energy implies that the Ali- 
Bodmer phenomenological aa potential cannot describe 
the ground state of 12 C with a compact shell- model like 
structure. 

On the other hand, the present 3a model interacting 
via the aa RGM kernel gives enough binding and a large 
overlap with the compact shell- model like component. 
Table III lists the results of such Faddeev calculations 
for the ground state of the 3a system with and with- 
out the Coulomb force. The aa RGM kernels are gener- 
ated from the VN1, VN2, and MN forces. When the 
Coulomb effect is included, the cut-off Coulomb force 
with Rc = 10 fm is employed. In the last column in 
Table III, C(o4) implies the overlap amplitude of the 3a 
bound-state function with the SU 3 (04) shell-model con- 
figuration. We find that all three effective NN forces 
yield binding energies comparable with the experimen- 
tal value \E™ P \ = 7.275 MeV, although the result of the 
MN force is a little too large. The dominant component 
of these 3a ground states is the SU3 (04) shell-model 
configuration. 

In Table IV we compare the 3a ground-state energies 
E3 a , predicted in the present three-cluster Faddeev for- 
malism, with those obtained by fully microscopic calcula- 
tions, Et^ 1 . We find that the present three-cluster equa- 
tion gives 3a energies which are only 1.5 - 1.8 MeV higher 
than those of the fully microscopic 3a RGM or GCM cal- 
culations. This implies that the three-cluster exchange 
effect, which is neglected in our three-cluster formalism, 
but is present in the fully microscopic three-cluster RGM 
kernel, is attractive in nature, and is not as large as the 
repulsive three-body force claimed necessary in the semi- 
microscopic 3a models [2, 4]. This is mainly because the 
3a model space used by these authors does not exclude 
the 3a Pauli-forbidden components accurately, unlike the 
one used in the present Faddeev formalism. 

In Tables III and IV, we also find that the three-range 
MN force gives a somewhat large overbinding of 2 - 4 
MeV, if the 3a energy E^ a is measured from the 3a 



TABLE IV: Comparison of the 3a ground-state energies, pre- 
dicted by the present model (E^a) and by fully microscopic 
calculations {E-J^ 1 ). The experimental value is E c £f — —7.275 
MeV. The present model is the Faddeev calculation using the 
aa RGM kernel, including the cut-off Coulomb force with 
Rc = 10 fm. The heading E^ implies the internal energy of 
the (0s) 4 a-cluster with the h.o. width parameter v, E to t the 
total energy from the RGM ([53] for MN and [38] for VN2) or 
GCM ([39] for VN1) calculations, and E^ 1 = E to t - 3E™\ 



Force 


v (fm" 2 ) 


77 int. 
£j a 


-Etot 


crfull 


E[ ia 


VN1 


0.2515 


-27.0 


-87.9 


-6.9 


-5.37 


VN2 


0.275 


-27.3 


-89.4 


-7.5 


-5.97 


MN 


0.257 


-23.9 


-83.0 


-11.4 


-9.59 



TABLE V: Kinetic- and potential-energy contributions to the 
three-a energy i?3 a , calculated from {Ho) = 2 (3e2 Q — E Sa ) 
and (V) = 3{E3 a — 2e2a)- The shell-model (04) component, 
C( 4), is large if {Ho) is large. 



Force 


£2a 


E?, a 


{Ho) 


{V) 


C(04) 


VN1 


10.778 


-5.366 


75.402 


-80.768 


0.875 


VN2 


9.503 


-5.969 


68.958 


-74.927 


0.790 


MN 


13.481 


-9.592 


100.068 


-109.660 


0.971 



threshold. The decomposition of the 3a energy to the 
kinetic-energy and potential-energy contributions in Ta- 
ble V implies that this overbinding is due to the large 
cancellation between these two contributions. In this re- 
spect, it is interesting to note that the a clusters with 
v = 0.257 fm~ 2 (which gives the correct rms radius 
r a = {2,/As/u) = 1.48 fm [58] for the simple (0s) 4 a- 
cluster) give less binding in the framework of the orthog- 
onality condition model (OCM) [8]. If the h.o. constant 
parameter v is small, a proper treatment of the aa ex- 
change kernel seems to be essential in order to obtain a 
large binding energy of the 3a ground state. This is rea- 
sonable since the large overlap of two a-clusters implies 
the importance of nucleon exchange effects. 

B. a a A Faddeev calculation 

For a detailed description of the aaA bound states in 
the Faddeev calculation, it is important to make sure 
that the result is converged with respect to the following 
three conditions: 

1) convergence with respect to the momentum dis- 
cretization points, 

2) convergence with respect to the extension of partial 
waves included, 

3) convergence with respect to the cut-off radius Rc, 
when the cut-off Coulomb force is included. 

Among them, the Coulomb effect is the most difficult, 
since the T-matrix of the full Coulomb force is divergent 
at the diagonal part and the strong oscillation in the 
momentum representation in the cut-off Coulomb case 
does not lead to the correct answer, unless the numer- 
ical angular-momentum projection of the aa Coulomb 
kernel (especially the direct Coulomb term) is accurately 
performed. As to the partial waves, we can easily enu- 
merate all possible angular-momentum states of ^Be for 
the 17 = 0+ ground state with J = 1/2 and the L n = 2+ 
excited state with J = 5/2 and 3/2 in the LS cou- 
pling scheme. If no Aa spin-orbit force is introduced, 
the J = 5/2 and 3/2 excited states are degenerate and 
the LS'-coupling scheme is more efficient than the jj- 
coupling scheme to reduce the number of channels cou- 
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pled in the calculation. In the following, the angular- 
momentum truncation is specified by A max -£ lmax values 
for the aa and Aa pairs. For example, D-P in the 
ground-state calculation implies a 4-channel calculation 
and G-G in the L v = 2 + calculation a 19-channel calcu- 
lation. The largest model space adopted is /-/, which is 
an 11-channel calculation for L" = + and a 28-channel 
calculation for L" = 2+ . Note that the variational calcu- 
lation in Ref. [7] uses a rather restricted model space; i.e. 
a three-channel calculation with A max = 2 and ^i max = 0, 
although the meaning of angular-momentum truncation 
is a little different from ours. For the momentum dis- 
cretization points, we find that the energy change due to 
the increase of ni-n 2 -n 3 is very much Rc dependent. It 
is usually positive if we go from m-n2-ri3=5-5-5 to ni-n 2 - 
n 3 =10-10-5 when the Coulomb force is not included, but 
it turns out negative when Rc = 10 fm and 14 fm. This 
implies that the Faddeev calculation without Coulomb 
force usually overestimates the binding energy, if the 
number of momentum discretization points is not large 
enough. Since the cut-off Coulomb kernels are oscillating, 
too small a number of momentum discretization points 
such as in ni-n 2 -n 3 =5-5-5 case is dangerous when Rc is 
very large like Rc = 10 fm and 14 fm. The orthogonal- 
ity to the Pauli-forbidden states also deteriorates when 
the number of momentum discretization points is too 
small. The squared norm of the Pauli-forbidden compo- 
nents contaminating the total wave function is typically 
10 -5 - 10~ 6 when m-n 2 -n 3 =5-5-5, but is improved to less 
than 10~ 13 for m-n 2 -ri3=10-10-5. In this paper, we will 
mainly show the results of ni-n 2 -n 3 =10-10-5, since the 
energy gain by further extension to ni-n 2 -n 3 =15-15-5 is 
usually less than 1 keV, when the cut-off Coulomb force 
with Rc = 10 <~ 14 fm is included. 

The energy gain of the ground state, AE, and that of 
the self-consistent e 2ct value by the increase of the maxi- 
mum angular-momentum values, A max -£i max , are shown 
in Table VI in the cases when we use the VN2 or MN 
forces for the aa interaction and the SB force for the Aa 
interaction. In these calculations the cut-off Coulomb 
force with Rc — 6 fm is employed. If the S'-wave calcula- 
tion is extended to include the D-wave, the energy gain is 
about 1 McV for VN2+SB and 1.2 MeV for MN+SB. The 
energy gain mainly comes from the partial-wave compo- 
nent with l\ = £2 = 1 of the a - ^He channel. The effect 
of the partial wave £\ = £2 = 2 is rather small; i.e., about 
50 (VN2) - 60 (MN) keV. Needless to say, the exact en- 
ergy gain largely depends on the character of the AN odd 
force. The ground-state energy is further improved by 7 
(VN2) - 5 (MN) keV and 0.03 (VN1) - 0.0 (MN) keV, 
according to the extension to the G- and /-wave, respec- 
tively. On the other hand, e 2a is improved by 160 (VN1) 
- 281 (MN) keV, 5 - 6 keV and 0.5 - 0.6 keV, according 
to the extension to the D-, G- and /-waves, respectively. 
In conclusion, partial waves up to the D-wave are suffi- 
cient within 10 keV accuracy. If we wish to have a 1 keV 
accuracy, we need to take into account at least up to the 
G-wave. This implies that the partial-wave truncation 



TABLE VI: Energy gain for the ground state (AE) and that of 
the self-consistent E2a value (As2a) in keV, for the extension 
of the maximum angular-momentum values, A m ax-^i m ax- The 
cut-off Coulomb force with Rc — 6 fm is included. 



Force VN2+SB 





AE (keV) 


Ae 2a 


(keV) 




5-5-5 


10-10-5 


5-5-5 


10-10-5 


S-S — > D-P 


-954 


-954 


165 


160 


D-P -> D-D 


-50 


-50 


5 


5 


D-D -> G-G 


-7 


-7 


6 


6 


G-G -> I-I 


-0.03 


-0.03 


0.6 


0.6 


Force 




MN+SB 






AE (keV) 


Ae 2a 


(keV) 


n\-ni-m 


5-5-5 


10-10-5 


5-5-5 


10-10-5 


S-S -> D-P 


-1165 


-1172 


287 


281 


D-P -> D-D 


-57 


-58 


7 


7 


D-D -> G-G 


-6 


-5 


7 


6 


G-G -> I-I 


-0.1 


-0.0 


0.5 


0.5 



in the Faddeev formalism is very efficient and the result 
converges very rapidly, according to the increase of the 
partial waves taken into account. 

Table VII shows the Rc dependence of the two-a en- 
ergy E( 8 Be), the self-consistently determined e 2a , the 
three-cluster ground-state energy E'(^Be), the A sepa- 
ration energy defined by B\(\Bc) = E( s Bc) + M\ — 
E(\Be), and the expectation value of the Aa Hamil- 
tonian, e\ a , when the momentum discretization points 
with ni-n 2 -n 3 =10-10-5 and the partial waves up to I-I 
are used in the MN plus SB model. The energy increase 
(and the accumulated one) due to the increase of Rc 
is also shown with the plus sign in the second (and the 
third) row. We find that the ground-state energy E(\Be) 
increases by 1.621 MeV when we move from Rc = to 
Rc = 6 fm, which is larger than 1.127 MeV calculated 
for the free two-a bound state. This seems to be nat- 
ural, since the two-a subsystem is more compact in the 
^ Be system. The energy increase in the self-consistently 
determined e 2a values is 1.435 MeV, which is about 200 
keV smaller than the energy increase in E(\Be), but is 
still larger than in the free two-a bound state by about 
300 keV. This observation is a good example that our 
self-consistent procedure of determining e 2a is reason- 
ably functioning. It is interesting to note that this large 
Coulomb effect in the three-body ground state; i.e., about 
1.4 times larger than in the two-a system, is characteris- 
tic for the increase of Rc from to 6 fm. For the range 
from Rc — 6 to 10 fm, just the opposite is true and the 
energy increase in the three-body ground state (85 keV) 
is smaller than in the two-a system (> 133 keV). This is 
apparently because the free aa relative wave function is 
more widely spread than the correlated aa relative wave 
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function in the A Be ground state. The tendency of £ 2 a 
falls just into the middle of these two extremes. By using 
this feature, we can easily estimate the full Coulomb ef- 
fect in the E( A Be) ground state. We find that the result 
with Rc = 10 fm is accurate within a 1 keV error both 
for E( A Bc) and e 2ct . From Table VII, the final result for 
the MN+SB potentials is 

£( A Bc) = 27.35 - 34.18 = -6.837 MeV , 

e 2a = 19.46 - 18.27 = 1.181 MeV , 

s Aa = 9.215 - 7.954 = 1.261 MeV , 

c (40 ) = 0.695 . (30) 

Here we have shown the kinetic-energy and potential- 
energy contributions separately in each energy, and C( 40 ) 
is the overlap amplitude of the A Be ground-state wave 
function with the shell-model (40) wave function. [Note 
that the sum of the £ 2ct potential energy and twice of the 
s\ a potential energy is the potential energy of E( A Be), 
but this is not true for the kinetic-energy terms.] We 
have also carried out the similar analysis in the VN2+SB 
model. The converged result of the VN2+SB forces, in- 
cluding the cut-off Coulomb force with Rc = 14 fm, is 
given by 

E( 9 A Bc) = 21.21 - 28.09 = -6.879 MeV , 

s 2a = 13.64 - 12.99 = 0.649 MeV , 

e Aa = 8.264 - 7.548 = 0.715 MeV , 

c (40 ) = 0.569 . (31) 

If we compare this result with Eq. (30) for the MN force, 
we find that the energy gain by the more attractive VN2 
force is only 42 keV. This result is rather surprising, if we 
consider that the VN2 force gives a two-a bound state 
with energy E^a = —245 keV. The Aa interaction by the 
SB force is also more attractive than in the MN force 
case due to the different choice of the h.o. width pa- 
rameter v. In other words, the ground state energy of 
^Be is not much affected by the poor aa and Aa in- 
teractions, as long as we find a well converged value by 
taking enough partial waves and a large number of mo- 
mentum discretization points. On the other hand, the 
E2 a and e\ a values for the MN force are larger than 
those for the VN2 force by almost 500 keV. This may 
be related to the difference of v values in the two cal- 
culations. The smaller v value, 0.257 fm~ 2 , in the MN 
force calculation means more extended a-clusters than 
in the VN2 calculation (y = 0.275 fm -2 ), which implies 
in turn that the relative wave functions in the 2a and 
Aa subsystems should be more compact in the MN case. 
This can be confirmed by comparing the kinetic-energy 
contributions in E( 9 A Bc), e 2a and eAa in Eqs. (30) and 
(31). For example, the kinetic-energy contribution in e 2a 
is 13.64 MeV in the VN2 case, while in the MN case it 
has a much larger value 19.46 MeV. The compactness 
of the aaA relative wave function in the MN case is 
also reflected in the fact that C(4 ) is larger in the MN 



TABLE VII: Cut-off radius (Tic) dependence of the Coulomb 
energies in the two-a bound state energy _E( 8 Be), the two- 
a expectation value e2a, the three-body bound state energy 
E(\Be), the A separation energy _Ba(aB 6 )! and the Aa ex- 
pectation value SAa- Calculations are carried out by using 
ni-7i2-n3=10-10-5 and the partial waves up to I-I. The three- 
range MN force and the SB force are used with v — 0.257 fm -2 
for the h.o. width parameter of the a-clusters. The energy in- 
crease (and the accumulated one) due to the increase of Rc 
is also shown with the plus sign in the second (and third) 
row. The experimental A separation energy is B A xp ( A Be) = 
6.71 ± 0.04 MeV. The suffix "ext" stands for extrapolation. 



Rc (fm) 





6 


10 


14 oo 


£( 8 Be) 


-1.260 


-0.133 


> 








+1.127 


> +0.133 










(> +1.260) 




£'2a 


-0.384 


1.051 


1.180 


1.181 (1.181) ext 






+1.435 


+0.129 


+0.001 








(+1.564) (+1.565) (+1.565)ext 


E(lBe) 


-8.543 


-6.922 


-6.837 


-6.837 (-6.837)ext 






+1.621 


+0.085 


+0.000 








(+1.706) (+1.706) (+1.706)ext 


B A (lBe) 


7.283 


6.789 


> 6.837 




£Aa 


1.390 


1.228 


1.260 


1.261 (+1.261) cxt 



case, even though the binding energy is smaller. Com- 
paring the result in Eq. (30) with the experimental value 
£: cx P(^Be) = -6.62+0.04 MeV, we can conclude that the 
MN+SB combination overbinds the A Be ground-state en- 
ergy by 220 keV. This is partly because our SB potential 
is of the pure Serber type (u = 1). If we choose u = 0.82 
for the SB force, the combination with the present MN 
force and v = 0.257 fm~ 2 yields E( 9 A Bc) = -6.621 MeV. 
In this case, the ^He bound-state energy is —3.105 MeV. 

We list the results of various AN effective forces used 
by Hiyama et al. in Table VIII, when they are used in 
combination with the MN force for the aa RGM kernel. 
The calculations are carried out with ni-n 2 -n 3 =10-10-5, 
Rc = 10 fm, and the partial waves up to the G-wave, to 
obtain the converged results with the accuracy of 1 - 2 
keV. 

Table IX lists aaA Faddeev calculations for the 2+ 
excited state, including the cut-off Coulomb force with 
Rc = 14 fm. The momentum discretization points with 
ni-n 2 -n 3 =10-10-5 are employed. When the partial waves 
are restricted to D-S or S-D, the 2 + -state energy is lo- 
cated above the a + ^He threshold with the threshold 
energy —3.12 MeV. The listing therefore starts from the 
7-channel calculation with D-P. We find that the result 
is almost converged with I-I and Rc = 14 fm, within the 
accuracy of 1 keV. The final result for the 2 + excited 
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TABLE VIII: aaA Faddeev calculations for the L n = 0+ 
ground state, including the cut-off Coulomb force with Rc = 
10 fm. The aa RGM kernel is generated from the three- 
range MN force with u = 0.94687 and v = 0.257 fm" 2 for the 
h.o. width parameter of the a-clusters. The G-matrix based 
effective AN forces in Ref. [7] are used for the Aa interaction, 
by slightly modifying the short-range repulsive part to fit the 
A separation energy _Ba(aH 6 ) = 3.120 MeV. Partial waves up 
to A max are included in aa-A channel and those up to £i max 
are included in the Aa-a channel. The heading _E( A Be) is the 
three-body ground-state energy of A Be in the aaA model, e2a 
the two-Q expectation value determined self-consistently, and 
EAa the Aa expectation value, and C( 40 ) is the overlap with 
the shell model (40) wave function. 



force 


Amax ^lmax 


E( 9 A Be) 




EAa 


C(40) 




S-S 


-5.580 


0.909 


1.136 


0.606 


NS 


D-P 


-6.681 


1.122 


1.250 


0.683 




D-D 


-6.736 


1.133 


1.255 


0.686 




G-G 


-6.743 


1.132 


1.257 


0.686 




S-S 


-5.734 


0.764 


0.774 


0.579 


ND 


D-P 


-7.375 


1.136 


0.838 


0.693 




D-D 


-7.478 


1.159 


0.842 


0.697 




G-G 


-7.483 


1.157 


0.843 


0.697 




S-S 


-5.682 


0.802 


0.882 


0.587 


NF 


D-P 


-6.839 


1.009 


0.942 


0.666 




D-D 


-6.901 


1.021 


0.944 


0.669 




G-G 


-6.906 


1.020 


0.944 


0.669 




S-S 


-5.620 


0.862 


1.030 


0.599 


JA 


D-P 


-6.622 


1.022 


1.112 


0.667 




D-D 


-6.672 


1.031 


1.114 


0.669 




G-G 


-6.677 


1.031 


1.115 


0.669 




S-S 


-5.566 


0.915 


1.154 


0.606 


JB 


D-P 


-6.431 


1.027 


1.253 


0.664 




D-D 


-6.469 


1.034 


1.255 


0.666 




G-G 


-6.475 


1.033 


1.256 


0.666 



state in the MN+SB model is 

E = 29.47 - 33.40 = -3.926 MeV , 
e 2 a = 21.55 - 17.54 = 4.013 MeV , 
e A „ = 9.481 - 7.930 = 1.551 MeV , 
c (40) = 0.645 . (32) 

If we compare Eqs. (30) and (32), we find that the 3 MeV 
excitation energy of the 2 + state mainly comes from an 
increase of the two-a kinetic energy (2 MeV) and from 
the two-a potential energy (1 MeV). This clearly shows 
the rotational nature of the ground + and excited 2+ 
states, composed of the two-a cluster structure with a 
weakly coupled A. 

Table X summarizes the present results with the MN 
force for the aa RGM kernel. The SB result shows the 
overbinding of the \Be ground-state energy by about 220 



TABLE IX: Same as Table VIII, but for the L n = 2+ excited 
state with Rc = 14 fm. 



force 


/imax-tlmax 


-C'IaPCJ 


S2a 




c (40) 




n p 
\ > i 


^ 7Q7 

— o. / y / 


^ Q87 


1 ^98 


U.U^IO 


SB 


D-D 


—3.874 


4.014 


1.536 


0.645 




P P 


— o.yzu 




1 ^o 


U.U4:0 




I-I 


— 3 Q9fi 
o. y^u 


inn 


1 

l.JUl 


U.U4:U 




n p 
\ > i 


— O.I uu 


^ Q90 


1 ^1 fi 
l.Olo 


u.uoy 


NTS 

IN O 


Pi Pi 


3 779 
— o. / 1 z 












O.OOl 


3 Q49 


1.544 


U.U11 




I-I 


O.OOl 


O.iJ'rkO 


1.544 






n p 
\ > i 


A 177 


A 097 


1.1 ou 


n fi/i8 


NT) 

IN 1 J 


Pi Pi 


A ^18 
<±.u lo 


A 071 


_L . iO'l 


U.UOl 




P-P 


*±.ooo 


A OfifS 


1 1 X7 

I.IOI 


U.UJ1 




I-I 


— A ^1 


A 0fi7 


_1_ . lOO 


U.UJ1 




n p 
\ ) i 


— o.oOo 


^ 89^ 


1 99^ 

l.ZZO 


U.UO 1 


NF 


D-D 


o. yoo 


o.oo ± 


1.226 


fl f\1Q 
u.uoy 




P P 


o.yoi 


3 RAQ 




fl fi^Q 
u.uoy 




T T 
1-1 


3 Q81 


3 RAQ 




fl fi^Q 
u.uoy 




D-P 


-3.645 


3.805 


1.380 


0.635 


JA 


D-D 


-3.710 


3.827 


1.385 


0.637 




G-G 


-3.762 


3.825 


1.401 


0.637 




I-I 


-3.762 


3.826 


1.402 


0.637 




D-P 


-3.460 


3.775 


1.507 


0.632 


JB 


D-D 


-3.510 


3.792 


1.511 


0.633 




G-G 


-3.568 


3.793 


1.535 


0.634 




I-I 


-3.568 


3.794 


1.535 


0.634 



keV and too small excitation energy of the 2 + excited 
state by about 130 keV. Table X also shows a compari- 
son with the results by Hiyama et al. [7] for the G-matrix 
based effective AN forces. We find that their results 
are a little lower than our results by about 70 - 90 keV. 
Since their calculation is a variational calculation using a 
smaller model space than ours, this is not a convergence 
problem of the variational calculation. A possible reason 
is the difference between OCM and RGM in the aa part. 
They used aa OCM, while ours is aa RGM. The OCM 
usually gives more attractive results than the RGM. In 
fact, it is well known that 3a OCM usually gives a larger 
binding energy than the 3a RGM for the ground state of 
the 3a system [59]. A small difference in the exchange 
term of the Aa folding potential may also contribute to 
this difference. 

If we arrange the effective AA~ forces in Table X in the 
order of more attractive nature, we find 

ND (-7.483) > NF (-6.906) > SB (-6.837) 
> NS (-6.742) > JA (-6.677) > JB (-6.474) . 

(33) 
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TABLE X: Summary of the ground-state energy £ gr (0+) and 
the 2 + excitation energy E x (2 + ) in MeV, calculated by solving 
the Faddeev equation for the aaA system in the LS coupling 
scheme. The aa RGM kernel is generated from the three- 
range MN force with u = 0.94687 and v = 0.257 fm" 2 for the 
h.o. width parameter of the a-clusters. 



Vajv 


J5 gr (0+) (MeV) 


£ x (2+) (MeV) 




present 


Ref. [7] 




SB 


-6.837 




2.911 


NS 


-6.742 


-6.81 


2.912 


ND 


-7.483 


-7.57 


2.930 


NF 


-6.906 


-7.00 


2.925 


JA 


-6.677 


-6.76 


2.915 


JB 


-6.474 


-6.55 


2.907 


Exp't 


-6.62 ±0.04 


3.024(3) 


[30, 31] 






3.067(3) 



The experimental value —6.62 ± 0.04 MeV is located be- 
tween JA and JB. However, this does not mean that the 
Jiilich potentials JA and JB are the most correct AjV 
interactions. It is well known that the spin-spin central 
terms of these Jiilich potentials are completely wrong and 
that they fail to reproduce the observed energy spectrum 
of the and ^He systems [60]. As for the 2 + excitation 
energy, all the results in Table X are between 2.90 - 2.93 
MeV. They are too small by 110 - 130 kcV with respect 
to the average value 3.04 MeV of the two resonances re- 
cently observed by 7-ray spectroscopy [30, 31]. Since the 
experimental error bars are at most ±40 keV even in the 
(K,ir) reaction [26], this is a meaningful disagreement. 
It would be interesting to examine the Is splitting of the 
5/2+ - 3/2 + states, by introducing a small AjV spin-orbit 
force predicted by our quark-model interaction. 

In order to show that the present aa RGM kernel gives 
a better result than simple aa potentials, we show in Ta- 
ble XI some results of aaA Faddeev calculations using the 
Ali-Bodmcr potential, ABd [12], and the Buck, Friedrich, 
and Wheatley potential, BFW [11]. In these cases, there 
needs no self-consistent procedure to determine e 2a - We 
only use the SB potential for the Aa interaction, since 
results with other effective AjV forces are easily evalu- 
ated from the above discussion in the case of the aa 
RGM kernel. In these a-particle models, we customarily 
use h 2 /M a = 10.4465 MeV • fm 2 and e 2 = 1.44 MeV • fm. 
The momentum discretization points with ni-n 2 -n 3 = 15- 
10-5 are employed. For the aa Coulomb potential, the 
folding potential of the cut-off Coulomb force with the 
(0s) 4 shell-model wave function is used with Rc = 10 fm. 
The h.o. width parameter of the (0s) 4 a-cluster for this 
Gaussian folding is v = 0.27127 fm" 2 in the ABd case 
and v = 0.257 fm" 2 in the BFW case. In the ABd case, 
this v value corresponds to the Coulomb-force parame- 
ter P = V3/(2 x 1.44) = 0.6014 fm" 1 and the a rms 



radius, r a = ('i/A^/v) = 1.44 fm. Since this v value 
is also used for the a-cluster folding for the AjV poten- 
tial, the Aa bound-state energy i?(^He) is a little shifted 
from the fitted experimental value —3.12 MeV. [The dif- 
ferent h 2 /Mjsr value also affects this difference.] Since 
the energy change is only about 0.06 MeV, we do not 
readjust the potential parameters of the AjV force. In 
the BF W ca se, the v value, 0.257 fm -2 , corresponds to 
(3 = v/4jV3 = 0.58538 fm" 1 and the rms radius of the 
a-cluster, r a = ^/3/(2(3) = 1.48 fm. In this case the dif- 
ference of the AjV bound-state energy, 0.054 MeV, from 
-3.12 MeV is solely from the different h 2 /M N value. The 
bound-state solutions of the BFW potential are used for 
the pairwise Pauli-forbiddcn states. The elimination of 
the Pauli-forbiddcn components from the three-body to- 
tal wave function is always inspected by calculating their 
squared norm, which is of the order of 10" 13 . 

We find that the acvA ground-state energy by the ABd 
potential is lower than the result of the MN force in 
Eq. (30) by 0.3 MeV. Note that even in this case the 
energy gain from the higher partial waves than the S 
wave is appreciable, i.e., 0.7 MeV. This implies that the 
S-wave assumption adopted by Filikhin and Gal [34] is 
not valid. They used a little different version of the Ali- 
Bodmer potential (type (a) with 125 MeV modified by 



TABLE XI: aaA Faddeev calculations for the L w = 0+ 
ground state by the Ali-Bodmer (ABd) [12] and Buck, 
Friedrich, and Wheatley (BFW) [11] aa potentials. The SB 
AjV force is used for the Aa interaction. The cut-off Coulomb 
force is included at the nucleon level with Rc = 10 fm. The 
h.o. width parameters of the a-clusters are assumed to be 
v = 0.27127 fm" 2 (ABd) and v = 0.257 fm" 2 (BFW). The pa- 
rameters h 2 /M a = 10.4465 MeV -far 2 and e 2 = 1.44 MeV -fin 
are used. Partial waves up to A max are included in the aa- 
A channel and those up to ^i max in the Aa-a channel. The 
momentum discretization points with ri\-n2-n->, =15-10-5 are 
employed. The Aa bound-state energy E( A He) for the SB 
AjV force is given in the first column. 



ABd±SB 





Amax ^lmax 


E(iBe) 


S2a 




c (40) 




S-S 


-6.409 


0.970 


-0.503 


0.466 




D-P 


-7.091 


1.013 


-0.532 


0.497 


-3.183 


D-D 


-7.147 


1.013 


-0.526 


0.499 




G-G 


-7.153 


1.018 


-0.518 


0.498 




I-I 


-7.153 


1.018 


-0.517 


0.498 


BFW±SB 


E(lUe) 


Amax-^lmax 


E(iBe) 




SAa 


C(40) 




S-S 


-5.544 


0.861 


1.776 


0.630 




D-P 


-6.971 


1.147 


1.973 


0.724 


-3.066 


D-D 


-7.038 


1.155 


1.979 


0.728 




G-G 


-7.043 


1.161 


1.979 


0.728 




I-I 


-7.043 


1.161 


1.980 


0.728 
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120 MeV) and obtained E( 9 A Be) = -6.55 MeV in the S- 
wave approximation. We expect an energy gain of about 
0.7 MeV from the higher partial waves and their result is 
overbound, in comparison with the experimental value, 
-6.62 ± 0.04 MeV. In Table XI, we find that the BFW 
potential gives a better result than the Ali-Bodmer force, 
but the energy is still lower than in the MN force case 
by 0.2 MeV. In this case we find that the effect of partial 
waves higher than the S wave is quite appreciable, i.e., 
— 1.5 MeV. This is of course due to the inner oscillation 
of the relative wave function between the two a-clusters 
in the aaA ground state. The shell-model like (40) com- 
ponents are about 0.7 in amplitude, which is appreciably 
larger than C(4 ) ~ 0.5 in the Ali-Bodmer case. 

IV. SUMMARY 

The three-cluster Faddeev formalism using two-cluster 
resonating-group method (RGM) kernels opens a way 
to solve few-baryon systems interacting via quark-model 
baryon-baryon interactions without spoiling essential fea- 
tures of the RGM kernel; i.e., the non-locality, the en- 
ergy dependence proportional to the exchange normaliza- 
tion kernel, and the existence of pairwise Pauli-forbiddcn 
states in some specific channels. In this paper, we have 
applied this formalism to three-cluster systems involving 
a-clusters; i.e., the 3a and aaA systems. These systems 
involve all of the above three features for the microscopic 
interactions between composite particles. In particular, 
the aa interaction is a prototype of composite-particle 
interactions, in which the fully microscopic RGM cal- 
culation is easy and very successful. It, however, in- 
volves a somewhat complex kernel structure composed of 
three non-trivial Pauli-forbidden states, and the energy- 
dependence of the interaction is rather strong in the 
Pauli-allowed model space. In the present Faddeev for- 
mulation, the Pauli-forbiddcn components between pair- 
wise clusters are completely eliminated from the total 
wave function of the three clusters. This can be achieved 
by introducing a special type of RGM T-matrix calcu- 
lated from the two-cluster RGM kernel, which satisfies 
the T-matrix version of the orthogonality conditions to 
the relative motion between two clusters. The on-shell 
and half off-shell properties of the T-matrix are just the 
same as those of the ordinary T-matrix. This RGM T- 
matrix involves a relative energy of two clusters as a pa- 
rameter, which is determined sclf-consistently by calcu- 
lating the expectation value of the two-cluster Hamilto- 
nian with respect to the total wave function resulting 
from the Faddeev equation. The Faddeev equation us- 
ing T-matrices is equivalent to the pairwise orthogonality 
condition model (OCM) of three-cluster systems, inter- 
acting via two-cluster RGM kernels. A nice point of this 
formalism is that the underlying nucleon-nucleon (NN) 
and hyperon-nucleon (YN) interactions are more directly 
related to the structure of three-cluster systems than in 
the models assuming simple two-cluster potentials. 



We have first applied the present formalism to the 
ground state of the 3a system by using three different 
types of effective NN forces, the two-range Volkov forces, 
No. 1 (VN1) and No. 2 (VN2), and the three-range Min- 
nesota (MN) force. The three-range MN force reproduces 
the S-, D- and G-wave aa phase shifts quite well in 
the simple (0s) 4 -model of the a clusters. The compar- 
ison with the 3a RGM calculation has shown that the 
present three-cluster formalism using only the aa RGM 
kernel gives a good approximation to the microscopic 3a 
model. The difference of the ground-state energies pre- 
dicted by these two models is less than 2 MeV. The effect 
of the antisymmctrization among three a-clusters, which 
is neglected in our formulation, is attractive and is not 
so large, as long as the Pauli-allowed model space of the 
3a system is properly treated. It is also shown that the 
three-range MN force gives a lower ground-state energy 
than the two-range VN1 and VN2 forces, resulting in 
a somewhat large overbinding of 2 - 4 MeV, if the 3a 
ground-state energy is measured from the 3a threshold. 

The application to the ^Be system has proved that 
our three-cluster formalism is soundly extended to the 
systems with two identical clusters, in addition to the 
systems of three identical clusters like the 3a system and 
the triton system. Here we have introduced a new ef- 
fective AiV force, called the SB force, which is made 
from the quark-model predictions of the AN phase shifts 
by using an inversion method based on supersymmetric 
quantum mechanics [42]. The SB force consists of two 
simple two-range Gaussian potentials which reproduce 
the low-energy behavior of the 1 S'o and 3 Si AN phase 
shifts predicted by AN-~£,N coupled-channel RGM cal- 
culations using the model fss2 [46]. Since any central 
and single-channel effective AN force leads to the well- 
known overbinding problem of ^He by about 2 MeV [50], 
the attractive part of the 3 Si AN potential is reduced by 
about 10 % to reproduce the empirical A-separation en- 
ergy, S^ xp (iHe) = 3.12 ± 0.02 MeV. The odd-state A7V 
force is assumed to be zero (pure Serber type). In addi- 
tion to this SB force, we have also used the effective AiV 
forces in Ref. [7] for comparison. The Aa interactions are 
generated from these AN effective forces by the folding 
procedure with respect to the (0s) 4 h.o. wave function of 
the a clusters. 

In the aaA Faddeev calculation, enough partial waves 
up to AMax = ^lMax = 6 arc included both in the aa and 
Aa pairs since the relative wave functions between two 
a-clusters are oscillating at least in the relative S- and 
D-waves. The detailed analysis shows that the partial 
waves up to the D-wave are sufficient if we do not mind 
a 10 keV inaccuracy. If we wish to obtain a 1 keV ac- 
curacy, we need to take into account at least up to the 
G-wave. This implies that the partial wave truncation is 
very efficient even in the present Faddeev formalism. The 
energy gain due to partial waves higher than the S-wave 
is about 1 MeV for the VN2 force and 1.2 MeV for the 
MN force, when these aa interactions are used in com- 
bination with the SB force for the Aa interaction. The 
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Coulomb effect between the two a-clusters is included 
by a cut-off Coulomb force at the nucleon level. The 
cut-off radius, Rq = 10-14 fm seems to be sufficient 
for a 1 - 2 keV accuracy. In the present formalism, the 
structure change of two a-clusters inside A Bc is clearly 
identified by calculating the kinetic-energy contribution 
in the two-a expectation value e 2a . The comparison of 
the Coulomb contributions in the aa bound state, ei a 
and the A Be ground state with respect to the change of 
Rc is very useful to measure the compactness of the two- 
a configurations in various environments. It is confirmed 
that the + ground state and the 2 + exited state of A Bc 
are well described by the contracted two-a cluster struc- 
ture with a weakly coupled A-particle in the dominant 
S-wave component. In the present calculation using only 
central forces, the three-range MN force and the SB po- 
tential with the pure-Serber character can reproduce the 
ground-state and excitation energies of A Be within an ac- 
curacy of 100 - 200 keV. The results in Ref. [7] based on 
the OCM framework are also confirmed within 100 keV 
accuracy. On the other hand, the simple a-particle model 
using the Ali-Bodmer aa potential, ABd [12], and the 
OCM using the deep Buck, Friedrich, and Wheatley aa 
potential, BFW [11], with bound-state Pauli- forbidden 
states give an overbinding of the A Be ground state by 
530 keV and 420 keV, respectively, when the SB force is 
used for the Aa interaction. Although these energies are 
rather similar, the effect of partial waves higher than the 
S-wave is very different; i.e., 0.7 MeV in the Ali-Bodmer 
case and 1.5 MeV in the BFW case. It is natural that 
the aa interactions which yield an oscillatory behavior 
of the aa relative wave functions, like our RGM kernel 
and the BFW potential, need more partial waves with a 
larger energy gain. 

There are still many problems left for the future stud- 
ies. First of all, the readjustment of the 3 S attractive 
part of the SB AN potential is unsatisfactory from the 
viewpoint of using the fundamental baryon-baryon inter- 
actions. The Brueckner rearrangement effect in A Hc is 
fairly large even for the rather stable a-cluster [50]. In 
this sense, there is still no consistent description of the 
s-shell and p-shell hypernuclei even at the level of us- 
ing effective baryon-baryon interactions. A microscopic 
description of the Aa interaction may need a more de- 
tailed analysis based on the G-matrix theory, for which 
the folding formula given in Appendix B is very useful. In 
order to describe the A Be excited states realistically, we 
need to introduce the Aa spin-orbit force and solve the 
Faddeev equation in the jj-coupling scheme. The recent 
7-ray spectroscopy experiment [30, 31] indicates a very 
small spin-orbit splitting for the possible 5/2+ and 3/2+ 
resonances. It is interesting to examine the LS compo- 
nents of the quark-model AN interaction, in which the 
antisymmetric LS interaction (LS^) is by about a fac- 
tor two larger than in the Nijmcgen models. We expect 
a large cancellation between the ordinary LS interaction 
and this LS^ interaction. An interesting application of 
the present Faddeev formalism and the Aa T-matrix de- 



rived in this study is to the recent Nagara event [61] for 
A 6 A He. For the AA interaction, we can use the coupled- 
channel AA-SAf-EE T-matrix of the quark-model inter- 
action, fss2. A preliminary result [41] shows that fss2 is at 
present the only model which can reproduce an appropri- 
ate strength of the AA interaction, AB™£ = 1.01 ± 0.20 
MeV, deduced from the Nagara event. In a separate pa- 
per [22], we will also report another application of the 
present three-cluster Faddeev formalism to the hyper- 
triton system, in which the quark-model NN and YN 
interactions are explicitly used in the ANN and T,NN 
coupled-channel Faddeev formalism. In this system, a 
complete Pauli-forbidden state at the quark level exists 
in the AN-T,N subsystem. 
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APPENDIX A: REARRANGEMENT FACTORS 
OF THREE-CLUSTER SYSTEMS WITH TWO 
IDENTICAL PARTICLES 

In this appendix, we give a brief comment on the defi- 
nition of the rearrangement factors in the Dirac notation 
for general three-body systems with two identical parti- 
cles or clusters. The incorporation of spin-isospin degrees 
of freedom is essential for further applications to the hy- 
pertriton system [22] and the AAa system [41]. When 
one uses the Dirac notation, it is important to fix a co- 
ordinate system of the representation. We choose the 
standard system of the Jacobi coordinates with 7 = 3, 
and introduce the Jacobi coordinates in the momentum 
space, p = p 3 and q = q 3 . The other Jacobi coordi- 
nates p 1 , q 1 , etc. are similarly defined. For an arbitrary 
function i[)(p, q; 123) in 7 = 3, the effect of the cyclic 
permutation P(i23) of the symmetric group 5*3 is 

^(123)^(^2,92; 312) = P(l23)V(Pl,9i;231) 

= ^(p 3 ,q 3 ;123) , (Al) 

where 123 in ip(p, q; 123) stands for the spin-isospin vari- 
ables. For the transposition P(i2), Eq. (1) yields 

^ ( i2)V>(P3, <z 3 ; 123) = iK-ps, <z 3 ; 213) , 

-P(i2)-0(Pi,<?i;231) = i/>(-p 2 ,q 2 ; 132) , 
^W(P 2 , 925 312) = ^(- Pl , 9l ; 321) . (A2) 

Note that the momentum suffix a in p tt , q a , and the 
sign of p a , etc., are uniquely specified by the sequence 
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of 123. For example, (p 2 ,q 2 \(3) m Eq. (4) actually im- 
plies (p 2l q 2 ] 312|/3). In the following, we always use an 
abbreviated notation, ip = ip(p 3 , q 3 ] 123), in the stan- 
dard coordinate system 7 = 3. The total wave function, 
^(<3f3) 93) 123), in Eq. (2) is then compactly expressed as 

tf = V - P(12)P(123) 2 V + P(123) V • (A3) 

If we write the Faddeev equation in terms of ip and (p, it 
reads 

V = G T (1 ± P (12) ) P (123) V , 

p = GoT[P ( i23)V±P (2 3 )V ] , (A4) 

with T = Ti2 and T = P(i23)T'i3P(i23) _1 ) where Ti 2 and 
Ti 3 are the two-body T-matrices in the three-body space. 



The definition of the rearrangement factors in the 
Dirac notation is based on the assumption 

{p 3 , <z 3 ; i23|p ( i23)V) = ^(123) V(p 3 > 93; 123 ) 

= ^2^2^12) = J dp'dq' S(p' -p 2 ) S(q' - q 2 ) 
xP$*<p(p',j;123), (A5) 

where the function <p(p 3 , q 3 ; 123) is in Eq. (A3) and 
^(123) °P era tes only on the spin-isospin variables of 
f{P ' 1 q'i 123). With this ip in the (3 = 2 channel in mind, 
the standard procedure of the partial wave decomposition 
gives the following definition for the first-type rearrange- 
ment factor g-yp^q, q'',x): 



(p,qn\P(i23) 2 \p',<i\ P)3-2 = ^ { dx 5<J> S<yP ,i 1 +2^ W o?? q';x) = [ dpdfidpdq' (7|p,9;i23) 



x 5 \ p + q' 



m 2 



m 2 + mi 



q)S[p'-q- 



m 3 



m.3 + mi 



q') (p',q'; 123|/3). 



(A6) 



Here, pi and p 2 are given in Eq. (5c) with a general mass factor, ( = (4m 3 /mi). With this mass modification, Eq. 
(6) is valid with a more complete reduced rearrangement factor 



M) A G^iL) (STT z;i \\P^f\\STT z] (3) 



(LS— coupling) 





A 




I 






*i 


h 




e 


s 2 


j 






s'2 


32 




L 


S 


J 




L 


S 


J 



G \u)lL) (STT^WP^WSTT^P) (ij-coupling) 



• (A7) 



Here the square bracket implies the unitary form of the 9j coefficients and the quantum numbers are specified by 

(A8) 



\ 1 ) = \[(\l)LS]JJ z -TT z ) \f3) = \[(t 1 £ 2 )LS}JJ z ;TT z ) (LS-coupling) 
| 7 > = \[(\ Sl )I(£s2)j]JJ z ;TT z \(3) = WM^s'^JJ^TT,) (ij-coupling) 



The angular-momentum factors ^v^w/j^) with Ai = <~ A, A'i = <~ t\ are given by 

(2A + l)!(2^i + l)! 



(2Ai)!(2A 2 )!(2Ai)!(2A 2 )! 



XU 



112 



xY,(^20m\fO}(X' 2 0£ 2 0\f'0)(kOX 1 0\f'0)(kOX' 1 0\fO) 

fj' 



f L Ai I I /' L Ai I I Ai /' L 

X X 2 i J \l x X' 2 l 2 J I Ai / k 



(A9) 



where A = V2A + 1 etc. and A 2 = A — Ai, A 2 = l\ — X' x . In the spin-isospin reduced matrix elements of Eq. (A7), the 

permutation operator P^S) docs not cnan g e tne total spin and isospin values, S and TT Z . 

The other types of rearrangement factors are obtained in a similar way. First, the symmetry of the matrix elements 
yields 



(P,q,(3\P(123)\p',<l'n)2-3 = (p',q',l\P(123) 2 \P,<l,(3}3- 



(A10) 



The rearrangement factor for the matrix clement (<p\P(23)) If) needs a little care, since the mass assignment of the 
three particles is made in the standard Jacobi coordinates 7 = 3. We first use P( 2 3) = P(i23)P(i2)P(i23)~ an d write 
the matrix element as 



( ( fi\P(23)\ ( fi) =X] / d P3 d <l3 V*(P 2 ,<?2;312) ^(-Pi,Qi;321) . 

123 •* 



(All) 
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The corresponding rearrangement factor in the Dirac notation is given by 

(p,q,f3\P {23) \p',q',(3'h-l = V / dpdqdp'dq' (ftp, q; 123) 6 ( p + q' + ^ q) s(p' + q+ " 2 q' 



x Pg } 123|/3') , (A12) 

from which the results in Eqs. (8) and (9) are easily obtained. Note that the rearrangement factor Eq. (A12) is 
symmetric with respect to the interchange between p, q, (3 and p', q', (3' , since mi = m 2 . 



APPENDIX B: A USEFUL FORMULA FOR THE 
Aa BORN KERNEL 

The general procedure to calculate Born kernels of the 
s-shell clusters, developed in Ref. [62], can also be used 
to calculate the Aa Born kernel 

V(q f , qi ) = {e i( lf r \V\e^ r ) 

5 

where <p a is the internal wave function of the a cluster, 
£a is the spin wave function of the A particle and vij is 
an effective AN interaction. The essential part of this 
method lies in the correct treatment of the cm. motion 
which is handled by the procedure given in Ref. [63] . This 
method makes it possible to deal with the most general 
form of the AN interaction with non-static effects like the 
G-matrix AA interaction. In this method, V(qf,q i ) in 
Eq. (Bl) is calculated from an integral form of the GCM 
kernel through 

5 
3=2 

= (2_y e ^(q f 2 +q< 2 ) J dadb e -*<lf a+ ^ b G(a,b), 

(B2) 

with 

G(a,b) = J dR (Ma)1> a (0)\ 

5 

xJ2 v v\MR + b)ip a (R)) ■ (B3) 

3=2 

Here, 7G = (4 + Qv, 7 = 4<>/(4 + Q with C = M A /M N , 
and iPa(R) and ip a (R) are the h.o. shell model wave func- 
tions of A and a, centered at R, with the width param- 
eters C,v and v, respectively. First we calculate spatial 
integrals for the spatial part u in v\j = liywy. These 
four integrals with j = 2 - 5 are all equal because of 
the antisymmetric property of the a cluster. We need 
to calculate spatial integrals for u = u{x\ — x 2 ) and 
u = u(xi — x 2 )P r , which we call the direct term and 
the exchange term, respectively. It is important to note 



that the space exchange operator, P r , operates only on 
the single-particle coordinates X\ and x 2 , and does not 
exchange the A and N masses. The procedure to inter- 
change these masses Ma and Mjv simultaneously like in 
Ref. [7] leads to an erroneous expression [see Eq. (A.l) 
of [7]], which is apparently wrong since the RGM kernel 
(S(r — a)£_A(j) a \ X^-=2 v ij\d( r ~ &)£a should not involve 
the mass dependence. The correct expression is the one 
in which one sets Mm — M\ in their Eq. (A.l) [see Eq. 
(B8) below]. The most general form of the two-body 
AA matrix elements for the translationally invariant u is 
parameterized as 

(PiPMpM = j^5(P-P')u(k',q';P) , (B4) 

with p = (p l - CP 2 )/(C + 1), p = Pi + P 2 ( also v'-, p ' 
for pi, p 2 ), and k' = p — p', q' = (p + p')/2. For the 
matrix element Eq. (B4), the spatial part of the GCM 
kernel in Eq. (B3) is calculated to be 

G spacc (a, b) 

= Wf^{Hf) 2 JdPdk'dq'u(k',q';P) 

x ex p/-l^±ip3 _ J_C + 1 ( ,2 1 ,2\ 
XCXP \ 6K+1 2u C V 4 fe J 

+i(a b) ■ (V + + il(a + b) ■ k j . 

(B5) 

If we use Eq. (B5) in Eq. (B2), we can perform the inte- 
grals over a and b and obtain two delta functions. Thus 
we can perform the integrals over k' and q' and obtain a 
compact formula 

3 

V sp&cc (q f , qi ) = e-^ fe2 (J^j 2 J dP e-^ p2 
where k — q i q i and q = + q,j) /2. 
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For a simple local Gaussian interaction, we find 




for u(r) = e- Kr P r . (B7) 



Then the P integral is carried out easily and we obtain 

2 

for u(r) = e~ Kr , 

v ^ qi) = \jVTY-J cxp l ^ 

-|^, 2 } for „(r)-e-^P P . (B8) 
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